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Abstract 

In this paper we study the relaxation singularities of human motor system at aging. Our purpose 
is to examine the structure of force output variability as a function of human aging in the time and 
frequency domains. For analysis of experimental data we have developed here the statistical theory 
of relaxation of force output fluctuation with taking into account the effects of two relaxation 
channels. The first of them contains the contribution of short-range correlation whereas other 
relaxation component reflects the effect of long-range correlation. The analysis of experimental 
data shows, that the general behavior of relaxation processes at human aging is determined by 
a complicated combination and nonlinear interactions two above stated relaxation processes as a 
whole. 
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I. INTRODUCTION 



Recent studies the aging of human neuromuscular system has focused on application of 
central theory regarding the age relatingwith the structure of behavioral and physiological 
variability due to loss of "complexity" . The term of "complexity" is connected to 

the broader concept including the fractals and dynamics in disease and aging. The concept 
consider behavioral and physiological changes of the system's due to aberrations in their 
time and frequency structure. However it is necessary to specify that less attention has 
been given to the neuromuscular changes associated with the time-dependent and dynamic 
peculiarities of force control till now. 

In this study, we examine the time and frequency structure of healthy adult humans 
to determine a relaxation singularities, arising with aging of human motor systems. It is 
well known that aging tends to induce a range of performance decrements in human motor 
system. A common finding is that the amount of motor variability increases in the healthy 
aging adult over a broad set of tasks 

Although consideration of the amount of variability is an important indicator of the aging 
neuromuscular system, the structure of motor variability also provides significant insight into 



system organization and human motor control j^l, 0, 10, 11]. The presence of short-time and 
long-range correlations in neuromuscular fluctuations in healthy people has implications 
for understanding and mathematical modelling of neuroautonomic regulation. Here we have 
applied correlation analysis to assess the effects of physiological aging on correlation behavior 
of human motor system. The purpose of this paper was to consider the effects of aging and 
task demand on the relaxation structure of the signals of force output variability from the 
point of view of modern statistical physics. We consider the force-time series data as a 
discrete stochastic process and apply the statistical theory and the information measure of 



memory of non-Markov random processes in complex systems [12|, 14, 15]. Therefore 
we can use the notions, based on the theory of chaos and information, in our analysis of 
complex systems. This allows to apply in our consideration representations and notions on 
statistical short-time and long-time memory, information measures of memory, relaxation 
and kinetic parameters from statistical theory jl6, 17, 1^ . 
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II. METHODS 



A. Participants 

A total of 29 participants were assigned to three different age groups: young group (N=l 0; 
range: 20 - 24 yrs; mean: 22 + 1 yrs; 5 females and 5 males), old group (N=9; range: 64 - 
69 yrs; mean: 67 + 2 yrs; 4 females and 5 males), and older-old group (N=10; range: 75 - 
90 yrs; mean: 82 + 5 yrs; 5 females and 5 males). All of the participants were right hand 
dominant. The participants were familiarized with the purpose of the experiment and all 
participants gave informed consent to all experimental procedures, which were approved by 
the local Institutional Review Board. 

The three age groups consisted of moderately active individuals. Persons who were highly 
active were excluded from the study. Also, elderly persons who were considered frail were 
excluded from the study. Twelve of the participants in the two elderly groups were taking 
medication for the treatment of high blood pressure. The distribution of persons taking 
medication for high blood pressure was 8 in the older-old group and 4 in the old group. None 
of the elderly participants reported having a neuromuscular or neuropsychiatric disease, nor 
did any of the participants have diabetes. Also, twelve of the elderly participants reported 
having arthritis and each were taking medication for the condition. The distribution of 
persons reporting arthritis of the hand was 7 in the older-old group and 5 in the old group. 
All participants remained on their normal medications during testing. 

B. Apparatus 

Participants were seated in a chair with their dominant forearm resting on a table (75 
cm in height). The participant's dominant hand was pronated and lay flat on the table 
with the digits of the hand comfortably extended. The setup constrained the wrist and 
the third, fourth and fifth phalanges from moving. The elbow position remained constant 
throughout the experimental session. Through abduction, the participant's lateral side of 
the index finger contacted the load cell (Entran ELFS-B3, New Jersey), 1.27 cm in diameter, 
which was fixated to the table. The load cell was located 36 cm from the participant's body 
midline. Analog output from the load cell was amplified through a Coulbourn Type A 
(Strain gauge Bridge) S72-25 amplifier at an excitation voltage of 10 V and a gain of 100. A 
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computer controlled 16-bit A/D converter sampled the force output at 100 Hz. The smallest 
increment of change in force the A/D board could detect was .0016 N. The force output was 
displayed on a video monitor located 48.6 cm from the participants' eyes and 100 cm from 
the floor. According to previous work from laboratory (D.A.V.), the display-to-control gain 
was set at 20 pixels per 1 N change in force for each participant. 

C. Procedures 

During the initial portion of the experiment, the participant's maximum voluntary con- 
traction (MVC) was estimated consistent with previous work \7\. Participants abducted 
their index finger against the load cell with maximal force for three consecutive 6 s trials. A 
60 s rest period was provided for each participant between each MVC trial. In each MVC 
trial, the mean of the greatest ten force samples was calculated. The means obtained from 
three trials were averaged to provide an estimate of each participant's MVC. Participants 
adjusted their level of force output to match a red target line (1 pixel thick) on the video 
monitor. Participants viewed online feedback of their performance in the form of a yel- 
low force-time trajectory that moved from left to right in time across the video monitor. 
They were instructed to match the yellow trajectory line to the red horizontal target line 
throughout each trial, and to minimize all deviations of the yellow line from the red line. 

Participants produced force at a constant force target. The constant target was a hori- 
zontal line displayed across the center of the video monitor. Participants produced force at 
5, 10, 20 and 40 % of their MVC under the constant force condition for two consecutive 25 s 
trials at each force level. A rest period of 100 s was provided between each force trial. The 
order of the force and target conditions was randomized across all participants. 

III. DATA ANALYSIS 

The force-time series data were conditioned by the following methods. First, force data 
were digitally filtered by using a first-order Butterworth filter with a low-pass cutoff fre- 
quency of 20 Hz. All data processing and subsequent time and frequency analysis were 
performed by using software written in Matlab. The analysis of force output concentrated 
on two primary problems: 1) calculation of the information measure of memory of force 
variability; 2) study of the correlation and relaxation structure of force time variability. 
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A. The information measure of memory for force variability. 



For comparison the relaxation time scales of the initial TCF a(t) and memory functions of 
the ith order Mj(t) (they will be described below in Eqn. (3)) we use here the first measure 
of memory that is the statistical non-Markovian parameter. Originally the non-Markovian 
parameter, characterizing the degree of non-Markovity of an arbitrary relaxation process, 
was introduced for analyzing the irreversible phenomena in condensed matters 16|, |17| • The 
relaxation times of initial TCF (the existence duration of correlations in considered system) 
and memory functions of the ith order (the duration of existence of memory) can be deter- 
mined as follows: r a = At J2f=o a (^j)> t m± = ^J2f=o Mi(tj). The simplest criterion for 
the quantitative estimation of non-Markovity in the given relaxation process is determined 
as: £1 = Tu/tm-l- When £j > 1, the relaxation time of the memory function of the first 
order is much smaller than the relaxation time of the initial TCF. In this case the process 
is characterized with a very short memory, in the limit e — > oo it is a Markovian process. 
Particularly, our work 18] proposes the way of the transformation from non-Markovian ki- 
netic equations to Markovian in case, when the non-Markovian parameter tends to infinity. 
Decreasing the parameter e\ determines the relative memory lengthening and strengthening 
of non-Markovian time effects. Thus, the presented quantitative criterion characterizes the 
degree of non-Markovity and strength of memory in the underlying relaxation process. 

Later 0] a conception of non-Markovian parameter spectrum e and markovization 
depth for nonequilibrium processes in disordered condensed matter was introduced. These 
parameters are related to fundamental properties of the system as well as the memory func- 
tion, the memory life time and demarkovization of the process by means of the initial TCF. 
The spectrum of non-Markovian parameter {e} = {e\, £2, ■ ■ ■ , £ n -i} is a set of dimensionless 
£i values: 



£l — Ta/TMi, £2 — Tm-l/tMv, ■ ■ ■ , £n-l — 7"Af„_i PM n ■ 

Here n is a relaxation time of memory function of ith order, the number % defines the 
relaxation level. 

The quantitative information measure of memory of force variability was assessed by 
calculating the statistical spectrum of non-Narkovity parameter (NMP). The equation: 
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where % = 1,2, 3, ... is the number of relaxation levels in kinetic description 12J of the time 
series, allow to calculate NMP Si from the power spectra of memory function of ith order 
Hiiuj). The mathematical procedure for a finding of ^{uj) consist in the calculation: 



Hi(w) 



n-1 



Y^Mitfr) cos(jujT) 

j=0 



(2) 



where Mi(t) (t = jr) is the memory function of z'th order, u = 2tt/t, t = 0, 01s is the dis- 
cretization time, n is the number of signals in time series. Memory functions, phase portraits 
of the dynamical orthogonal variables, and set of dynamical orthogonal variables of junior 
orders were calculated by the methods of well-known statistical theory of stochastic dis- 
crete non-Markov processes in complex systems with applications to cardiology, seismology, 



physiology etc. The full details of the theory have been described elsewhere 
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B. The correlation and relaxation analysis of force variability. 

The structure of correlation, relaxation and memory processes in force variability will 
be examined by using time and frequency analysis of correlation and studied relaxation 
processes. Correlation and relaxation rates for the short-range Rf {Rf = |Aj|,i = 1,2) 
and long-range Rf\R^ — |Aj \ 1 ^ 2 , % = 1,2) relaxation (correlation) were used to consider the 
force signal. We shall calculate a general relaxation time in a frame of relaxation theory 
developed below. 

Let's consider the structure of the initial time correlation function a(t) of the force output 
time series signals. It is convenient to consider the chain of interconnected equations for the 
TCF a(t) in a frame of well-known Zwanzig'-Mori's kinetic theory 12, I3I. [bll. flo| : 



' da(t) 
dt 
dMi(t) 
dt 



Xxaty-Ax f a{t -T)Mi(r)dT, 
Jo 



A 2 Mi(t) - A 2 / M 1 (t - r)M 2 (r)dT, 



(3) 



Here, the parameter |Aj| and lAjl 1 / 2 define the local and memory channel relaxation (corre- 
lation) rates on the ith relaxation levels, correspondingly. Parameters Aj has the dimension 
of relaxation rate, and Aj has a dimension of a square of frequency. 

The structure of correlation and relaxation processes in force variability will be examined 
by using time and frequency analysis of correlations and studied relaxation processes. Cor- 
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relation and relaxation rates for the short-range (R^ = | Ai|, i = 1,2) and Ion g-range 
Rf ) (Rf ) = \A i \ 1 / 2 ^i = 1,2) relaxation (correlation) were used to consider the force signal. We 
shall calculate a general relaxation time in a frame of relaxation theory developed below. 
Let's consider a partial solution of system (3) when we can use the time-scale invariance 



idea in nonequilibrium statisticalw physics of condensed matter 24 



25, MM 



M 2 {t) = M 1 {t) (4) 

on the second relaxation level. It means the approximate equality of memory life-times on 
the first and second relaxation levels. The similar condition is carried out in many concrete 
cases. The condition (4) apply on our situation as we shall see later. Using the Laplace 
transform on the normalized TCF a(t): 

poo poo 

a(s) = / dte~ st a(t), M,(s) = / dte'^M^t), (5) 
Jo Jo 

memory functions Mi(t) and M 2 (t), one can obtain solution of Eqns. ([3]): 

M x {s) = ^{-(s + R2) ± v / (^ + J R 2 ) 2 + 4A 2 }, (6) 

a(s) = {s + R 1 + 2^[-(a + #2) + 
v / Gs + i? 2 ) 2 + 4A 2 ]}- 1 , 

where Ri = = |Aj|, i=l, 2. If we determine the general relaxation time tr by the 
relation tr = 3? lim s ^o 5(s), where 3? means real part, the resulting general equations for tr 
consequently can be written in the following way: 



Ai 



r R = {Ri + ^[\/ 4A 2 + R 2 - 3a]}" 1 . (7) 



if 4A? + Ro > 0. Then we receive the formula: 



rn = Rl " 2 , (8) 

(i?i-^R 2 ) 2 + ^|4A 2 + i?!| 

for A 2 < 0, 4A 2 + R 2 2 < 0. 

The equations ([8]) are our new theoretical results and they allow to calculate a 
general relaxation time tr for the various relaxation scenario. One can see from Eqns. (JTj), 
(jSJ) that general relaxation behavior is defined by complicated nonlinear interactions and 
combinations of short-range and long-range relaxation processes on the first and second 
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relaxation levels. Therefore, the general behavior of relaxation time tr is determined by 
a complicated combination of relaxation rates R^\ R^\ IA1J 1 / 2 and IA2I 1 / 2 on these two 
interconnected relaxation levels. Calculating these five relaxation parameters : tr, t± = 
Ri 1 , = -R2" 1 , Ti = IAxI" 1 / 2 and = | -A. 2 1 — 1//2 we can receive a rather detailed and 
specific singularities of relaxation processes in a complex systems. 



IV. RESULTS 



Figure 1 shows force output from young (B), old (C), and older-old (D) participants at 



20% MVC at the constant (tl) and sine wave (t2) targets. The amount of 



brce variability 



was examined by calculating the root mean square error (see, for details 28j|). 

Figure 2 depicts the values of short-range relaxation parameter on a first relaxation level 
Ai (Ri = I Ai| is the corresponding relaxation rate) with following name structure. 
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The file naming includes the experimental design AGE GROUP (3) x FORCE LEVEL (4) 
x TRIALS (2): 

si, s2, s3 (group number): sl-young (20-24 years), s2-old (64-69 years), s3-oldest old (75-90 
years); 

cl, c2, c3, c4 (force level): cl-5%, c2-10 %, c3-20 %, c4-40 %; 
tl, t2 (trial number): tl-trial 1, t2-trial 2. 

Each file contains 25 s force data (in Newtons) collected at 100 Hz {r = 0, 01s). A more 
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detailed description of the experiment and data collection is found in 

Figure 2 shows numerical values of parameter Ai for all three groups si, s2 and s3. In 
Figure 2 the box lines at the lower quartile, median and upper quartile values. The whiskers 
are lines extending from end of the box to show the extent of the rest of the data. Outliers 
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Figure 1: Force output from young (B), old (C), and older-old (D) participants at 20% MVC at 
the constant (tl) and sine wave (t2) targets. 

are data with values beyond the ends of the whiskers. We should add still, that first eight 
boxes correspond to si group, next eight boxes correspond to the second (s2) group and 
last ones correspond to the third (s3) group. One can calculate from Fig. 2 that the 
mean value of parameter Ai for first (si) group is Ai = —0,0203 r" 1 , for second (s2) group 
Ai = —0,0110 r" 1 , and for third (s3) group of participants Ai = —0,00763 t^ 1 . Similar 
distinction in rates and times of the short-range relaxation is quite explainable with the 
physiological points of view. We shall note, that the greatest distinction in young and old 
groups achieves 3,47 times. From physical point of view it means a more fast relaxation of 
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Figure 2: Relaxation parameter due to short-range correlations Ai (or short-range relaxation rate 

(s) 1 

R\°> = | Ai|) for the all studied groups at first relaxation level in units of r , where r = 0, 01 s 
is a discretization time. The details of description and processing of statistical data see in text 
of the paper. Data testify for fast relaxation for young (R[ s) = 0,02 t- 1 ), more slow relaxation 
for old (R[ s) = 0,01 r" 1 ) and oldest old (r[ s) = 0,007 r 1 ) groups. In a whole a difference for 
short-range relaxation rate for ratio of young / old constitute 1,8 times, for young / oldest old is 
2,66 times and for old / oldest old is 1,44 times. In a separate cases last ratio constitute 3,47 and 
more times! 

the force fluctuation for young and more slower relaxation for old and oldest old. 

Figure 3 show a similar behavior of the second short-range relaxation parameter A2 for all 
studied groups of participants. Authentic appreciable distinction (in 2,11 times) in relaxation 
rates for young for force levels cl, c2 and c3, c4 pays on itself attention. Average values 
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Figure 3: Relaxation parameter A2 due to short-range correlation (rate of short-range relaxation 

(s) 1 

R\ = IA2I) for the all studied groups at second relaxation level in units if t . For young group 
one can see a sharp decreasing of R 2 an d A2 at hight force levels (c3, c4). In the old group it is 
seen a steady decreasing of relaxation rate at raising of force level as for constant well as well for 

(s) 

sine wave force output task. In oldest old group steady decreasing of R 2 * s observed for constant 
force output task whereas for sine wave force output task a decreasing is definite less. 

i?2^ in units of r" 1 in young group (0,778) and old (0,795) almost coincide with each other. 
Distinctions between force levels cl, c2 and c3, c4 in group oldest old is almost twofold. The 
greatest relaxation rate (R 2 = 1, 056 r _1 ) is registered in young group at low force levels 
(cl, c2). The least relaxation rate is observed as for young group (R { 2 s) = 0,50 r" 1 ) at high 
force levels (c3, c4), and for oldest old group s3 {R^ = 0,36 r _1 ) at high force levels (c3, 
c4). 
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Figure 4: Relaxation parameter A3 (rate of short-range relaxation R 3 = | A3 1 ) due short-range 
correlations for all studied groups at third relaxation level in units of r" 1 . One can observe almost 
invariable behavior of A3 for a si, s2 and s3 groups. For a young group with enhancement of force 
level minimum is observed for force level c3 (20 %). 

The behavior of the third relaxation parameter A3 is displayed on Figure 4. Relaxation 

(s) 

rate R 3 = | A3 1 on the third relaxation a level is appeared, approximately, identical for the 
all three age groups. 

Relaxation parameters A x and A 2 , connected with long-range time correlations, are shown 
in Figures 5 and 6 in units of r -2 . Relaxation time (r x = jAil 1 / 2 ) for the first relaxation 
levels are almost identical to all three age groups (si, s2 and s3). It testifies that long-range 
correlations mechanism of force relaxation, practically does not depend on age (average t[ 
for groups si, s2 and s3, are equal, respectively, 0,98 r; 1,16 r and 0,98 r). Relaxation 
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Figure 5: Relaxation parameter due to long-range correlation Ai (appropriate long-range relaxation 
rate R± = IA1J 1 / 2 ) for si, s2, s3 groups at first relaxation level in units of r~ 2 . As parameter Ai 
has a dimension of a square of frequency, parameter R^p = |Aj| 1//2 has a dimension of relaxation 
rate in units of r -1 , where r = 0,01 s is discretization time. One can see a steady increasing of 
this rate in four times for young group at constant force output task and twofold increasing for sine 
wave output task. For old group R^' is steady for various force levels, in fact. Similar behavior of 
R^p is remained for oldest old group in process of increase of force levels. 

parameter A 2 turn out as negative for all three age groups si, s2 and s3. This testify to 
a change of relaxation mode at the transition from the first on second relaxation levels. 
Relaxation times r!p for groups si, s2 and s3 changes slightly (they are equal 0,96 r ; 0,79 
r and 1,11 r, correspondingly). It can testify upon the weak age changes in long-range 
relaxation mechanism of force output fluctuations. 
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Figure 6: Relaxation parameter due long-range correlation A2 (of long-range relaxation rate = 
(A) 1 / 2 ) for all studied groups si, s2 and s3 at second relaxation level in units of r . One of the 
specific peculiarities of this case is a negative sign of parameter A2. It means a change of relaxation 
mode on the second relaxation level. For all studied groups si, s2 and s3 we have approximately 
similar relaxation behavior. Consequently from the physical point of view one can to suppose that 
contribution to relaxation due to long-range correlations is almost steady with aging. 

In Fig. 7 one can see a values of non-Markovity parameter £i(0) of the first relaxation 
level on zero frequency for three age groups si, s2 and s3. Appreciable interest represent, 
that in young group (average value E\ = 51) and oldest old group (si = 55, 5) the chaoticity 
and randomness is approximately the same. Simultaneously it is obvious, that the greatest 
chaoticity {e\ = 107) appears in group old participants s2. Fig. 8 presents statistical non- 
Markovity parameter on the second level £2(0) for the all age groups. One can note, that 
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6,(0) 




Figure 7: First point of the statistical non-Markovity parameter £i(0) on zero frequency u = 
as the information measure of memory for groups si, s2 and s3. Values of £i(0) is approximately 
same for young {e\ ~ 51) and oldest old groups and E\ ~ 107 for old group. It means a similarity 
of quantitative measure of memory and randomness (existence of Markov effects) for these two 
groups (young and oldest old) and increase of this measure for group s2. 

everywhere £2(0) is almost identical and equal unity. From results of Fig. 8 it follows, that a 
condition of applicability of time-scale invariance idea and approximate equality of memory 
life-time on the first and second relaxation levels (see, Eqn. (2)) is carried out with the 
high degree of accuracy. Because of the experimental data submitted on Fig. 8, one can use 
Eqns. (6) - (7) for calculation of general relaxation time tr. 

Fig. 9 presents a values of time tr calculated according our theory, Eqns. (6), (7), for all 
studied groups of participants. Mean relaxation time for young group (si) at enhancement 
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Figure 8: Behavior of the second point of non-Markovity parameter £2(0) in all studied groups si, 
s2 and s3. One can see identical values of £2(0) ~ 1 with fine precision. Namely similar behavior 
let us to use a time scale-invariance idea and apply a solution (6)-(7) of kinetic equations (l)-(3) 
for analysis of experimental data. This figure constitute an experimental basis of our solution of 
kinetic equations for TCF and MF. 

of force level (from cl, c2 to c3, c4) increase from value 42,5 r up to value 56,07 r. For the 
group of young mean m is 49,3 r, where as for the group of old and oldest old mean tr 
is equal 91,1 r and 99,8 r, correspondingly. It means, that a general relaxation time tr at 
ageing increases on two times, approximately. 

Fig. 10 shows the phase portraits in a planes of four junior orthogonal dynamic variables 



Wo, Wi, W2 and W3, calculated according theory 



3, HQ 



for experimental file named 



slclt2 (B4) as an example. These portraits appears symmetrical in planes: (Wi,W2; 
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Figure 9: General relaxation time tr, calculated in accordance with Eqns. (6), (7) of our theory for 
all studied groups in units of r. Aging appears as an approximate doubling of tr and slowdowningof 
a whole relaxation process including as a short-range and well as a long-range correlations. 

(Wi,W3) and (W2, W3). Slight deviation from this symmetry one can see for planes of 
junior variables (Wo, Wj), i=l ,2,3. Phase clouds are dense and concentrated. Fig. 11 
present phase portraits for file name slc4tl of young participant (B4) at highest force level 
(40 %). One can note a significant multifold swelling of volume of all phase clouds. This 
testify the noticeable increasing of chaoticity of motor force activity for this case. 

Figs. 12 and 13 depict a change of the structure of phase clouds for old participant (C2) 
for force levels clt2 and c4t4, as an example, correspondingly. From Fig. 12 one can find 
the effect of condensation of all phase clouds. It corresponds to lowest value of the first 
non-Markovity parameter (ei(0) ~ 18) and weakly marked chaoticity of force fluctuation. 
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Figure 10: Phase portraits in the planes (Wi,Wj) i,j = 1,2,3 of junior orthogonal dynamical 
variables Wi for participant B4 of young group at force level cl and sine wave force output t2. One 
can characterize these portraits by a dense nucleus and small dissipation of phase points on the 
planes. Similar behavior is typical for healthy young human. 

At the transition to other force level (c4, 40 %) from Fig. 13 one can see the remarkable 
increasing of the volumes of all phase clouds and their appreciable asymmetry. It lead to 
increasing of chaoticity of force fluctuation almost 7 times (^i(0) = 120). 

Phase clouds for oldest old participants (D8) for force levels c3tl and c4t2 are shown 
on Fig. 14 and Fig. 15, correspondingly. One can observe a change of distribution of 
phase points from consolidated ones (for planes (Wi,VT 2 ), (W^i,^) and (^2,^3)) to a 
more scattered distribution at transition from force state c3tl to state c4t2. Simultaneously, 
a slight amplification of asymmetry of phase clouds in planes of junior variables (Wo, Wj), 
j=l,2,3 occurs. Thus, a quantitative measure of memory £i(0) on the first relaxation level 
does not varies almost (it equal 24 and 27, correspondingly). Therefore, the degree of 



18 




0.15 
0.1 
0.05 


5 -0.05 
-0.1 
-0.15 



d) 




-0.2 1 ' 1 

-0.3 -0.2 -0.1 
W, 



0.15 

0.1 e) 

0.05 ■ 

c<r o ■ 

.It*. 
CO 

5 -0.05 
-0.1 ■ 
-0.15 



-0.2 ' ' ^ 

0.1 0.2 -0.3 -0.2 -0.1 

W, 




0.15 
0.1 
0.05 

<sr o 

CO 

5 -0.05 
-0.1 
-0.15 



0.1 0.2 



-0.2^ 
-0.2 




-0.1 



0.1 



0.2 



W 2 [x- 



Figure 11: Enhancement of force level (c4) for participant B4 from young group reduce the scat- 
tering and appearance asymmetry of distribution of phase points. Therefore, raising of force level 
and force capacity lead to the swelling of phase clouds and to increase of the disordering effects in 
phase space. 

Markovization and chaoticity for these two cases does not change. 

Figs. 16,17 and 18 demonstrate a frequency dependence of the first non-Markovity pa- 
rameter E\(uj) for young (B5) old (C2) and oldest old (D2) participants in all force levels as 
an examples. From Fig. 16 for participant D5 one can see that Markov effect at ultralow 
frequency remain constant at force level tl (constant). At sine wave (e,f,g,h) force output 
task one can observe a progressive and steady decline of £i(0) (from value of 50 on value 
of 8). It is connected with a reduction of Markov effects and steady amplification of non- 
Markov effects. Similar behavior of E\ (uj) reflects the appearance of slight robustness in force 
fluctuation for young participants in process of increasing of force level. A data for old (C2) 
presents opposite example, as can one see from Fig. 17. We see a drastic amplification of 
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Figure 12: Enhancement of force level (c4 instead of cl) for participant C2 from old group reduce 
the scattering and appearance asymmetry of distribution of phase points. Therefore, raising of 
force level and force capacity lead to the swelling of phase clouds and to increase of disordering 
dynamics in phase space. 

Markov effects in process of enhancement of force level for constant (a,b,c,d) and sine wave 
(e,f,g,h) force output task. This means a steady Markovization of a process, that is, a steady 
amplification of Markov random effects. From Fig. 18 one can note analogous behavior of 
£\{u) for oldest old (D2) at the constant (a,b,c,d) and sine wave (e-h) force output task. It 
testifies about a gradual and more distinctive display of non-Markov effects and robustness 
in force behavior for participant (D2). 
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Figure 13: Enhancement of force level (c4 instead of cl) for participant C2 of old group reduce to 
twofold swelling of phase clouds, remarkable stratification and scattering of phase points in phase 
space. 



V. DISCUSSION AND CONCLUSION 



The purpose of this study was to examine the structure in the time and frequency domains 
of force output variability as a function of human aging. 

For the analysis of force fluctuation we have used the statistical relaxation singularities 
of motor system related to the short-range and long-range time correlation. 

It is necessary to emphasize some fundamental points, following from the study. Force 
output dynamics for all three groups of healthy (young, old and oldest old) can be charac- 
terized rather a high level of chaoticity and randomness on the first relaxation level of force 
output fluctuation. For all files studied we have found out, that £\{u) changes in a wide in- 
terval of values (7-300). It allows to state a strong chaotization and sharp expressed Markov 
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Figure 14: Phase portraits in planes (Wi, Wj) i,j = 1, 2, 3 of junior dynamical orthogonal variables 
Wi for participant D8 of oldest old group at hight level force c3 (20 %) and constant force output 
task (tl). Dense phase clouds (d,e,f) are accompanied by the slight stratification (see, figs a,b,c) 
of phase clouds. 

effects for all healthy groups. As is known, beginning of illnesses, for example Parkinson 



disease (see, for instance, |23|) results to suppression of Markov effects and appearance of 
non-Markovity effects and robustness. 

The same level of a randomness observed for young {e\ ~ 51) and oldest old (e\ ~ 56) 
turned out to be unexpected. Thus, on average the rather highest level of a randomness (e\ ~ 
107) are registered in old group. These values of NMR (ei) signify healthy state for studied 
participants. Therefore it is possible to state with assurance that the information measure 
of mamory is the greater for old, but it is smaller for young and oldest old groups. Steady 
non-Markovity [e-i ~ 51) on the second relaxation level for all three age groups confirms our 
supposition for the advantage of time-scale invariance idea of relaxation processes of force 
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Figure 15: Twofold enhancement of force level (c4, 40 % instead of c3, 20 %) for participant D8 of 
oldest old group at sine wave force task reduce on minor change of the structure of phase portraits 
and phase clouds without any stratification and scattering of phase points in a phase space. 

output fluctuation on these level. 

The relaxation singularities of force output fluctuations consist in the following. At the 
first relaxation level a remarkable distinction of relaxation rates of force output fluctuation 
is observed. It means that relaxation related to the short-range correlation is a more fast at 
low force level. With the enhancement of force levels the relaxation is rather decelerated. 
Noticeably more slow relaxation appears in two old groups s2 and s3. Hence we can conclude 
that aging becomes apparent as the notably slowing-down of relaxation processes, connected 
with the effect off the short-range correlation. In the second relaxation level contribution of 
short-range correlation for all age groups are the same. At that we notice stably (more than 
2 times) decreasing of the relaxation rate with increasing of force level (from cl,c2 to c3,c4). 

The following fact attract our steadfast attention. One can note, that contribution of 
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Figure 16: Frequency dependence of the first point of non-Markovity parameter £\{uS) for partici- 
pant B5 of young group for 8 trials: at four force levels (cl, c2, c3, c4) and at two form (tl, t2) of 
force output task. Data demonstrate almost invariable behavior of £i(w) at constant force output 
task, and steady decreasing of quantitative measure of memory at sine wave force output task (for 
example, £i(0) decrease from value of 40 to value of 8). 

long-range correlations in relaxation rates is, approximately, identical for the all three age 
groups and it does not depend on age on the first relaxation level. However, in young group 
here specific features are observed. They are that here relaxations rates increase, on the 
average, at 1,5 times at hight force levels (c3 and c4). On the second relaxation level the 
similar picture is kept. 

For analysis of experimental data we have developed here the statistical theory of relax- 
ation of force output fluctuation with taking into account: first two relaxation levels and 
effects of two relaxation channels. One of the relaxation channel contains the contribu- 
tion of short-range correlation whereas other component of relaxation reflects the effect of 
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Figure 17: Frequency dependence of the first point of non-Markovity parameter £i(w) for partici- 
pant C2 of old group for 8 trials: for force levels (5 %, 10 %, 20 % and 40 %) and two form (tl, 
t2) of force output task. Values of quantitative measure of memory (ei(0)) increase almost at ten 
times at constant force output task, whereas it increase almost at five times at sine wave force 
output task. 

long-range correlation. The analysis of experimental data shows, that one can determine 
the general behavior of relaxation processes as a whole by a complicated combination and 
nonlinear interaction of these two above stated relaxation processes. 

VI. ACKNOWLEDGMENTS 

This work was partially supported by the RFBR (Grant no. 05-02-16639-a), Grant of 
Federal Agency of Education of Ministry of Education and Science of Russian Federation 
(Grant no. RNP 2.1.1.741). This work has been supported in part (P.H.) by the German 



25 




co [2 jc / t] 



co [2 tc / x] 




250 



200 



150 




100 



10"' 
co [2 7t / x] 



10" 
co [2 7t / x] 



Figure 18: Frequency dependence (a,b,c,d) of the first point of non-Markovity parameter £\(uj) 
for participant D2 of the oldest old group is related with the steady decreasing of quantitative 
measure of memory and randomness at enhancement of force level (5 %, 10 %, 20 % and 40 %) 
and constant force output task (tl). Similar behavior is observed for sine wave force output task 
(t2) also. Simultaneously, one can see low frequency range (from zero up to W 2 Hz) with the big 
values of £i(u;), that testify the existence a steady frequency region of Markov random effects. 
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